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We examine a class of operations for topological quantum computation based on fusing and measuring topo¬ 
logical charges for systems with SU( 2)4 or k = 4 Jones-Kauffman anyons. We show that such operations 
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I. INTRODUCTION 

The topological approach to quantum computation pro¬ 
poses the strategy of achieving fault-tolerance by utilizing 
nonlocal state spaces of topological phases of matter, which 
naturally protect encoded information from local perturba¬ 
tions 11 71. In particular, non-Abelian quasiparticles in such a 
system collectively possess a nonlocal state space which may 
be used to comprise topologically protected qubits. Topo¬ 
logically protected computational gates acting on such states 
can be generated through braiding operations obtained by 
transporting the quasiparticles around each other. Fusion of 
quasiparticles and topological charge measurement provide 
a mechanism for initialization and read-out of the computa¬ 
tional state, at which point topological protection is explicitly 
violated in order to extract the information from the topologi¬ 
cal state space. In Refs. QSSS , some additional operations 
utilizing fusion, measurements, and ancillary anyons were 
proposed for discrete gauge theories. It was later realized that, 
by utilizing ancillary anyons, one could generate all braiding 
transformations (also in a topologically protected manner) us¬ 
ing pair-wise topological charge measurements or by 

adiabatically tuning pair-wise interactions between quasipar¬ 
ticles lfl2h . 

Recently, the idea of utilizing fusion operations together 
with topological charge measurements in order to generate 
topologically protected computational operations has been 
further explored in Ref. [!13]. In that work, it was demon¬ 
strated that such operations could provide useful operations 
for Ising anyons/Majorana zero modes. These operations 
augment the computational power of the braiding gates for 
such quasiparticles, though they are still unable to yield a 
fully topologically protected computationally universal gate 
set. More generally, the utility of fusion and measurement op¬ 
erations for topological quantum information processing has 
been largely unexplored. 

In this paper, we advance the study of fusion and measure¬ 
ment operations in topological quantum computation by ex¬ 
amining SU( 2)4 and k = 4 Jones-Kauffman (JK 4 ) anyons, 
which are closely related theories. Potential physical realiza¬ 
tions of such anyons may occur as quasiparticles in fractional 
quantum Hall systems, such as k = 4 Read-Rezayi states III 
or second level Hermanns hierarchy states [Tjj], which, re¬ 


spectively, have filling v = 2/3 and 3/5. Such states provide 
possible candidates for physically observed Hall plateaus in 
the second Landau level ltl6l4l9ll . 

It is known that braiding transformations of SU( 2)4 or JK 4 
anyons are not computationally universal H. We will show 
that fusion and measurement operations augment the compu¬ 
tational power of braiding operations for such anyons, allow¬ 
ing us to produce a computationally universal topologically 
protected gate set. This is accomplished, in part, through 
the generation of an exact, topologically protected irrational 
phase gate. These results provide a compelling demonstration 
of the utility and value of studying fusion and measurement 
operations for the purposes of topological quantum computa¬ 
tion. 

We note that another approach to achieving topologically 
protected computational universality for SU( 2)4 anyons us¬ 
ing measurements and non-standard encodings of qutrits was 
considered in Ref. |f20f| . 

This paper is structured as follows. In Sec. [II] we provide 
a brief review of anyon models, topological charge measure¬ 
ments, and some details of the SU (2)4 and JK 4 anyon models, 
upon which the rest of the paper is focused. In Sec.HnJ we dis¬ 
cuss a number of different ways of encoding quantum infor¬ 
mation in collections of SU( 2)4 or JK 4 anyonic quasiparticles, 
and provide protocols for changing between different encod¬ 
ings using fusion and measurement. In Sec. lIVl we present the 
single qubit computational gates obtained from braiding oper¬ 
ations. In Sec.[V] we provide a simple fusion-based protocol 
that generates the NOT gate. In Sec. lVIl we provide a protocol 
for “topological qubit fusion,” which uses fusion and mea¬ 
surements to reduce the number of qubits and collapses the 
encoded state in a specific nontrivial manner. We explain how 
this protocol may be used to convert certain ancillary states 
into computational gates. In Sec. lVIIl we provide a protocol 
for generating a topologically protected irrational phase gate. 
This is carried out by first utilizing the operations (based on 
fusion and measurement) developed in previous sections to 
generate a specific ancillary state, and then utilizing topologi¬ 
cal qubit fusion to convert the ancillary state into the irrational 
phase gate. Combined with braiding gates, this provides a 
computationally universal set of topologically protected sin¬ 
gle qubit gates. In Sec. IVIII1 we provide a protocol for gen¬ 
erating a topologically protected controlled-Z gate, utilizing 
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the previously developed operations. In Sec. IIXI we make a 
few concluding remarks. In AppendixlAl we give general ex¬ 
pressions for the basic data of the SU(2)fc and JK^ anyons 
models, and tabulate the values of /•’-symbols and //-symbols 
of JK 4 anyons that are used in this paper. In Appendix [Bj we 
derive the probability of a symmetric ID random walk com¬ 
pletely avoiding negative-valued positions during the first n 
steps, which is useful for analyzing the implementation of the 
irrational phase gate. 

II. ANYON MODELS 


Most anyon models of interest have no fusion multiplicities, 
i.e. N° b = 0 or 1 , in which case the vertex labels /1 are triv¬ 
ial and can be left implicit. We only consider anyon models 
of this type in this paper, so we will drop these labels from 
now on. General states and operators are described using fu¬ 
sion/splitting trees constructed by connecting lines with the 
same topological charge. Charge lines corresponding to quasi¬ 
particles have definite values, while charge lines further down 
the fusion tree, corresponding to nonlocal degrees of freedom, 
permit superpositions of topological charge values. 

In this manner, inner products are formed diagrammatically 
by stacking the corresponding vertices, which gives 


In this section, we briefly review the basic fusion and braid¬ 
ing properties of quasiparticles (point-like localized excita¬ 
tions) in (2+l)D topological phases of matter, as described by 
anyon models (a.k.a. unitary braided tensor categories). For 
additional details, see Refs. Il2ll l22n and references therein. 

An anyon model has a finite set C of topological charges, 
which obey a commutative, associative fusion algebra 

a Xl) = Y^ N ab c : (1) 

cec 

where N° b are positive integers indicating the number of dis¬ 
tinct ways charges a and b can be combined to produce charge 
c. There is a unique “vacuum” charge, denoted 0, which has 
trivial fusion (and braiding) with all other charges (for exam¬ 
ple, JV ^ 0 = 5 ac ) and which defines the unique conjugate a of 
each topological charge a via N® b = Sab- 

A quasiparticle in the physical system carries a definite 
value of topological charge, because it is a localized object. 
The fusion rules indicate that overall topological charge of a 
collection of quasiparticles may take superpositions of differ¬ 
ent values, as long as the quasiparticles are well-separated and 
there is more than one fusion channel (i.e. N£ b and N£ b are 
nonzero for c c'). This gives rise to degenerate nonlocal 
state spaces associated with non-Abelian anyonic quasiparti¬ 
cles, which are topologically protected. 

More specifically, each fusion product has an associated 
vector space V£ b with dim V£ b = N£ b , and its dual (splitting) 
space V° h . The states in these fusion and splitting spaces are 
assigned to trivalent vertices with the corresponding topologi¬ 
cal charges. We write orthonormal basis vectors of the fusion 
and splitting spaces as 

(2) 

(STY-*.*«* **• ® 

where /j, = 1,..., N^ b . The normalization factors in terms of 
d a are included so that the diagrams will be in the isotopy- 
invariant convention, meaning bending lines and rotating sec¬ 
tions of diagrams only change amplitudes by unitary trans¬ 
formations (we refer the reader to Refs. i2lll22fl for details). 



This relation can be applied inside more complicated dia¬ 
grams. Note that this diagrammatically encodes charge con¬ 
servation. Since this formalism describes the states associated 
with anyonic quasiparticles (in a topological phase of matter), 
we require the inner product to be positive definite, i.e. d a are 
required to be real and positive. 

Associativity of fusion is represented in the state space by 
the /’-symbols, which (similar to (jj-symbols) provide a uni¬ 
tary isomorphism relating states written in different bases dis¬ 
tinguished by the order of fusion. Diagrammatically, this is 
represented as 

a b c a b c 

x/ -Em,, \x (5) 

d * d 

If the diagram on either side of this equation is prohibited by 
the fusion rules, the corresponding /’-symbol is set to 0 . 

The quantum dimension of topological charge a 

d a =d- a = \[Fr] 00 r ( 6 ) 

is also equal to the largest eigenvalue of the matrix N a defined 
by [N a ] bc = N£ b , and so describes how the dimensionality of 
the state space grows asymptotically as one introduces more 
quasiparticles of charge a (i.e. dim[l/ a -°] ~ </” when the 
number n of charge a quasiparticles is large). We also define 
the total quantum dimension of an anyon model to be V = 

The counterclockwise and clockwise braiding exchange op¬ 
erator of topological charges a and b are, respectively, repre¬ 
sented diagrammatically as 



The action of the braiding operator on the state space can be 
described in terms of //-symbols, which represent the unitary 
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operator for exchanging two anyons in a specific fusion chan¬ 
nel, and are obtained by applying the exchange operator to the 
corresponding trivalent vertices 



An anyon model is defined entirely by its N£ b , /-’-symbols, 
and //-symbols. The N° b must provide an associative and 
commutative algebra. The f'-symbols and //-symbols are 
constrained by the “coherence conditions” (also known as the 
“polynomial equations”), which ensure that any two series of 
F and/or R transformations are equivalent if they start in the 
same state space and end in the same state space Q Physi¬ 
cally, these consistency conditions are interpreted as enforcing 
locality in fusion and braiding processes. 

An important invariant quantity derived from braiding is the 
topological twist of charge a 



which is a root of unity. Another important invariant quantity 
is the topological //-matrix 


Sab = V - 1 




( 10 ) 


When S is unitary, braiding is non-degenerate and the theory 
is called a “modular” theory (MTC). 

The projector onto collective topological charge a of n 
anyons of definite charges a \,..., a n is given by 



where we sum over all possible fusion channels (with the same 
values in the bra as in the ket). 

For a modular theory (with unitary S'-matrix), we can write 
the projector of n anyons onto definite collective topological 
charge a by enclosing the charge lines of these anyons with an 
UJ a loop 



( 12 ) 


where the w a -loop is defined as 


co, = X>aS 0 * 


(13) 


A. Measurements 

There are two classes of topological charge measurements: 
local and interferometric. 

Local measurements are capable of measuring the topolog¬ 
ical charge ascribed locally to a single quasiparticle (or sim¬ 
ilar object) and can also be used to measure the topological 
charge of the fusion channel of a pair of quasiparticles. Such 
measurements are physically performed by measuring an ex¬ 
ternally measurable local quantity, such as energy or fermion 
parity, which is correlated with the topological charge value of 
a quasiparticle or small region. In order to measure pairwise 
fusion channels using such local measurements, one should 
interpret the diagram of a pairwise projector more physically. 
Specifically, to perform such a measurement, one must bring 
the pair of quasiparticles close to each other (or modify the 
system in some fashion that induces the same effect), essen¬ 
tially fusing them, so that one can apply a local measurement 
to the fusion outcome of the pair, before finally moving the 
quasiparticles apart to their original positions (splitting them). 

Interferometric measurements of topological charge 11221 
[241] are physically performed by sending probe quasiparti¬ 
cles through a device which generates coherent superpositions 
of the probe quasiparticles traveling through distinct paths 
around some region. The interference between different paths 
is correlated with the total collective topological charge con¬ 
tained in the encircled region, and this is detectable by mea¬ 
suring the probe quasiparticles exiting the interferometric de¬ 
vice. Such measurements are capable of distinguishing be¬ 
tween the topological charges ascribed to the overall fusion 
channel of the entire collection of quasiparticles contained 
inside the interference loop region of the device. These are 
non-demolitional measurements of topological charge, in the 
sense that, as long as the quasiparticles inside the interferome¬ 
ter are kept sufficiently well-separated (distances much larger 
than the correlation length), only their overall fusion chan¬ 
nel is collapsed by such a measurement. State information 
encoded higher in the fusion tree than the overall fusion chan¬ 
nel may be unaffected by the interferometric measurement. 
We emphasize that this effect cannot be achieved using local 
topological charge measurements, since bringing more than 
two quasiparticles close to each other will cause these fusion 
channels to evolve and decohere in some non-universal man¬ 
ner. 

The asymptotic operation of a generic ally tuned anyonic in¬ 
terferometer using probes that can distinguish all topological 
charge types will: (1) project the anyonic state onto the sub¬ 
space where the collective topological charge of everything 
contained inside the interferometry region (encircled by the 
probe paths) takes the definite value a, and (2) decohere all 
anyonic entanglement between the interior and exterior re- 
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gions of the interferometer | [25|| . The decohering effect (2) 
requires the use of density matrix formalism to describe, and 
it can be described using a superoperator acting on the den¬ 
sity matrix (as described in Refs. 11 111 l22l 12411 ) or, equiva¬ 
lently, by applying u>o loops to the density matrix in the di¬ 
agrammatic representation (as described in Ref. lEhll ). It is, 
however, possible to remove the decohering effect (2) with 
a protocol that essentially reconnects the severed topological 
charge line, which is projected onto charge a by the measure¬ 
ment, connecting the interior and exterior regions of the inter¬ 
ferometer III 3 Ll27h . In this manner, we can use interferometry 
measurements to generate only the topological charge projec¬ 
tion (1), which can equivalently be represented by inserting an 
bj a loop around the anyonic charge lines of the quasiparticles 
in the interior region of the interferometer. 

Thus, we can write the effect of both types of topologi¬ 
cal charge measurements using the standard von Neumann 
projective measurement formalism applied to state vectors. 
Specifically, a topological charge measurement for a state 'I') 
with measurement outcome equal to a will occur with proba¬ 
bility 


Pa = (^IrialT) = ||n a |T')|| 2 , (14) 


and transforms to the (normalized) post-measurement state 


I*) ^ 


n 0 |tt) 

l|n„|*>ir 


(15) 


Throughout this paper, we will always represent measure¬ 
ment projectors using cu-loops, regardless of whether the mea¬ 
surement used is local or interferometric (though when the u> 
loop encircles the charge lines of more than two quasiparti¬ 
cles, it necessarily requires an interferometric measurement 
to achieve). 


B. SU(2)4 and JK4 Anyons 


The set of topological charges for SU(2)4 and JK 4 anyon 
models is EH 


C = {0,1,2,3,4}. 
The fusion rules are given by 


0x0 = 0, 0x1 = 1, 

0x3 = 3, 0x4 = 4, 

1x2 = 1 +3, 1x3 = 2 

2x2 = 0 + 2 + 4, 2x3 = 1 

3x3 = 0 +2, 3x4 = 1, 


(16) 


0x2 = 2, 

1 x 1 = 0 + 2, 

4, 1x4 = 3, (17) 

3, 2x4 = 2, 

4x4 = 0 


(notice that a = a for all a ) and the quantum dimensions 


For JK 4 the topological twist factors are 

9q = O4 = 1 , 9\ = —0 3 = e l 4 , 82 = e 1 3 . ( 20 ) 

There are two other anyon models possible with these same 
fusion rules, and they are simply the complex conjugates 
SU( 2)4 and JK 4 of the these two theories (the only difference 
in the above data being that the twist factors are all complex 
conjugated). 

There are too many ^-symbols and i?-symbols to write 
them all out explicitly, but it is straightforward to compute 
them using the general expressions given in Appendices I A II 
and IA 21 We tabulate the F-symbols and /+symbols for JK 4 
that we will use for calculations in this paper in Anncndix IA 31 

In the following, we focus only on JK 4 , since these are com¬ 
pletely isotopy invariant (all bending and raising and lowering 
of lines at vertices are done freely) 1431 , which simplifies the 
calculations. In particular, this allows us to apply f-moves 
and /Amoves in any orientation within a diagram. The basic 
results that we derive for JK 4 will carry over to SU ( 2)4 anyons 
with only minor modifications. The fact that the results carry 
over follow from the fact that the two theories can be related 
by simply gluing a semion onto the odd integer charges. El 


III. ENCODING QUANTUM INFORMATION 

There are various ways in which one may encode quantum 
information in the nonlocal topological state space of non- 
Abelian anyons. Even when there is a single nontrivial non- 
Abelian topological charge in the theory, one may choose to 
encode states more densely or sparsely in the state space. 

The “standard encoding” of a topological qubit is given 
by four non-Abelian quasiparticles of the same topological 
charge q whose collective fusion channel (of the four quasi¬ 
particles) is fixed to the vacuum charge 0 , when a pair of 
charge q quasiparticles has two possible fusion channels. 
(More generally, if any two of the topological charges q have 
d possible fusion channels, then this would provide a qudit.) 


A. Encoding qubits in JK4 and SU(2)4 

For the JK 4 and SU( 2)4 theories, one may use q = 1 to 
form a standard encoding qubit, whose basis states are the two 
fusion channels a = 0 and 2. These are represented diagram- 
matically as 


1111 



a 


do — t/4 — 1, d\ — do — a/ 3 , d 2 — 2. ( 18 ) 

For SU( 2)4 the topological twist factors are 

0o = f ?4 = l, d 1 = -6 3 = e i T, 0 2 = e^. (19) 


We refer to this as the “111 1” encoding. 

Similarly, the standard encoding using q = 3 also provides 
a qubit with fusion channels a = 0 and 2. In contrast, the 
standard encoding using q = 2 provides a qutrit, whose basis 
states are the three fusion channels a = 0, 2, and 4. 
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Another qubit encoding, which will form the primary com¬ 
putational basis utilized in this paper, is given by four quasi¬ 
particles, two of which carry topological charge 1 and two of 
which carry topological charge 2, with the collective fusion 
channel of the four quasiparticles being fixed to vacuum 0. 
The topological charges of the two fusion channels encoding 
the qubit basis states depend on which pair of quasiparticles 
we fuse. If we fuse a 1 and a 2 quasiparticle, the two fusion 
channels are a = 1 and 3. We refer to this last qubit encoding 
as the “1221” encoding, and represent the basis states as 


12 2 1 



a 


If we fuse the 1 and 1 quasiparticles with each other, or the 2 
and 2 quasiparticles with each other, the two fusion channels 
are a = 0 and 2. (A pair of charge 2 quasiparticles may poten¬ 
tially have fusion channels 0, 2, and 4, however, the collective 
fusion channel restriction of this encoding excludes the charge 
4 fusion channel.) This encoding can be related to the 1221 
encoding via the /•’-move 

1221 1221 

= E K 22 L 

£>= 0,2 

so we do not view it as a distinct encoding, but rather as a 
change of basis. 

One may also consider more dense encodings of qubits. In 
particular, it may be useful to consider two qubits that are en¬ 
coded in the collective state space of six quasiparticles. For 
six quasiparticles with respective topological charges 1, 2, 2, 
2, 2, and 1, whose combined fusion tree is given by 


1 2 2 2 2 1 



a f, 


we have basis states for two qubits given by the possible topo¬ 
logical charge values a = 1,3 and b = 1,3. We refer to this 
as the “two qubit 122221” encoding. 

B. Switching from the 1111 encoding to the 1221 encoding 

In order to switch between different encodings of quantum 
information, one generally needs to utilize operations that in¬ 
volve fusion and measurement. (For a theory in which braid¬ 
ing is dense, one may use braiding operations to change how 
information is encoded. However, this approach may not be as 
efficient as the fusion and measurement based approach con¬ 
sidered here. Moreover, we are particularly interested in the¬ 
ories for which braiding is not dense.) 

We now describe a process that non-deterministically takes 
one from the 1111 encoding to the 1221 encoding. Starting 
from a qubit in the 1111 encoding, we follow the steps: 


1. Introduce an ancillary pair of charge 2 quasiparticles 
that are pair-produced from vacuum. 

2. Fuse one of the charge 1 quasiparticles with one of the 
charge 2 quasiparticles and measure the resulting fusion 
outcome topological charge x £ {1, 3}. 

3. Fuse the (new) charge x quasiparticle with one of the 
original charge 1 quasiparticles and measure the result¬ 
ing fusion outcome topological charge y £ {0, 2,4}. 

The above steps may be applied in many different fashions, 
e.g. where the ancillary quasiparticles are produced, which 
quasiparticles are fused, etc. Depending on the specifics of 
this process, one may need to spatially rearrange the quasi¬ 
particles, e.g. via braiding operations, during or after these 
steps in order to obtain the desired physical processes and final 
configuration. These details and choices involved in enacting 
this protocol (following steps 1-3) are most easily specified 
diagrammatically. For the measurement outcomes x and y, 
we consider the implementation of this protocol specified di¬ 
agrammatically by: 


IV 2 1 



In this diagram, we have added dotted horizontal lines as a 
visual aid for partitioning the diagram into three sections cor¬ 
responding to the three steps of the protocol. 

In step 2, the probability of obtaining either fusion outcome 
x = 1 or 3 will be p x = i. In step 3, the probability of 
obtaining a particular value of y will, in general, depend on 
the initial state of the qubit. Thus, the operation in step 3 will 
read out some state information and collapse the qubit state 
(at least partially, unless the initial state is in a basis state |a)). 
As usual with measurement processes in quantum mechanics, 
the state is re-normalized to unit norm. 

In order to obtain a state in the 1221 encoding, we require 
the final fusion outcome to be y = 2. When y ^ 2, the quasi¬ 
particles no longer support a degenerate state space, so they 
can no longer encode a qubit. In this case, we must discard or 
recycle the quasiparticles. As such, we should only use this 
protocol on ancillary qubits. 

Evaluating the diagram for y = 2, we obtain the transfor¬ 
mation in terms of qubit basis states of the initial 1111 encod¬ 
ing and the final 1221 encoding after fusion and measurement 
with outcomes x = 1 or 3 and y = 2 to be 

till 1221 

- E p tJ 

6=1,3 

a 5 

where 

P l: ) = ^[ F f 1 } lb [ F b 11X }a2 ■ ( 2 7 ) 
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(Here, we wrote the relation without the normalization con- This result is obtained from the following diagrammatic eval- 
stants on the diagrams, and will similarly drop overall con- uation 144] 
stants throughout the paper, whenever they are not important.) 

I 




EPI,. 





12 2 1 


b 


(28) 


We reemphasize that transformations pi x ^ are not unitary 
and may (partially) collapse the initial state. Hence, the 
overall factors here are not important, because one must re¬ 
normalize the post-measurement state anyway. Evaluating the 
^-symbols gives 


pi * 1 2 ) 


1 

7! 


1 V2 

1 4 

pi 3) _ J_ 

f) 3y/2 

U 4 

n 3V2 

u 4 

v/3 

1 \/2 
_ 4 


(29) 


where the columns are assigned values b = 1,3 and the rows 
are assigned values a = 0, 2. 

Thus, this implementation of the change of encoding proto¬ 
col transforms the initial qubit state 


|T)=^T a |a) (30) 

a—0,2 

in the 1111 encoding into the final qubit state 

p(P> l^r) 


Pi*) IT) 

P (x) |*> = E P ba* a \b) (32) 

a—0,2 
6=1,3 


T (31) 

2 



in the 1221 encoding. 

Notice that the outcomes x = 1 and 3 are related by appli¬ 
cation of a NOT gate 


X = 


0 1 

1 0 ’ 


(33) 


i.e. P^ = XP( 3 \ Thus, it is not important which measure¬ 
ment outcome was obtained for x during the protocol, as long 


as we are able to apply a NOT gate when needed to obtain the 
desired outcome. In Sec.|V] we will describe how to obtain a 
(deterministic) X gate using a simple fusion process. 

We observe that the diagram on the right hand side of the 
first line of Eq. (1281) simply displays the diagrammatic rep¬ 
resentation of the protocol in a more symmetric form. As 
such, it emphasizes that different possible protocols (i.e. fus¬ 
ing and measuring in different orders) may be used to produce 
the same resulting encoding changing operation. We will not 
go into detail on these, but should keep in mind that alternative 
realizations may potentially be advantageous. 


C. Switching between two qubits in 1221 encodings and a two 
qubit 122221 encoding 

We can switch the encoding of two qubits from two 1221 
qubit encodings to a two qubit 122221 encoding, with no state 
collapse, by using a “forced measurement” procedure 
m, as follows. Starting from two 1221 qubits (enumerating 
the quasiparticles 1-8 from left to right), we follow the steps: 

1. Measure the collective topological charge x £ {0, 2} of 
quasiparticles 4 and 5. If x ^ 0, go to step 2. If x = 0, 
go to step 3. 

2. Measure the collective topological charge y £ {0, 2} of 
quasiparticles 5-8. (The outcome is not important.) Go 
to step 1. 

3. Remove the (now ancillary) quasiparticles 4 and 5 in 
a manner that does not create entanglement, e.g. fuse 
them together into vacuum or transport them as a pair 
along the same path away from the other quasiparticles. 
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Steps 1 and 2 constitute a forced measurement, i.e. a “re¬ 
peat until success” measurement process, which may take n 
attempts to achieve the desired x = 0 measurement outcome. 
More specifically, the nth iteration of the measurement in step 
1 is given by 



where y± = 0 , and the nth iteration of the measurement in 
step 2 is 


12211221 



= K 1 ] 


XnVn+l 



so failure is exponentially suppressed in the number of at¬ 
tempts. In the notation of Refs. li, we write the forced 
measurement procedure (achieving success at the nth attempt) 
as 


n^ 45) = ni 4 n 5 l 0 n^ 678 )... n^ 6 J 0 8) . (38) 

We emphasize that the measurements made in steps 1 and 2 
are independent of the encoded state, as reflected by the inde¬ 
pendence of the measurement outcome probabilities on a or 
b. This indicates that these measurements do not collapse or 
otherwise alter the encoded state. 

Combining steps 1 and 2 with step 3, which just removes 
the now ancillary pair of quasiparticles, we obtain the change 
of encoding 


12211221 



a b 

1 2 2 2 2 1 


(39) 


a b 

which does not collapse or otherwise alter the two qubit state 
encoded in the fusion channels a and b. 

It should be clear that the inverse of this change of encod¬ 
ing, taking us from the two qubit 122221 encoding (enumerat¬ 
ing the quasiparticles 1-3, 6-8 from left to right, with positions 
4 and 5 left vacant for the introduction of ancillary quasipar¬ 
ticles) back to the two 1221 qubits, may be obtained through 
the steps: 

1. Introduce an ancillary pair of charge 1 quasiparticles 
that are pair produced from the vacuum, taking the va¬ 
cant positions of quasiparticles 4 and 5. 

2. Measure the collective topological charge y £ {0, 2} of 
quasiparticles 5-8. (The outcome is not important.) If 
y ^ 0, go to step 3. If y = 0, stop. 

3. Measure the collective topological charge x £ {0, 2} of 
quasiparticles 4 and 5. Go to step 2. 

This utilizes the forced measurement 



The probability of measurement outcome x n = 0 at the nth 
attempt is 


Po 




(36) 


^(5873) = n (5678) n (45) _ _ _ n (5678) ^ _ (40) 

IV. SINGLE QUBIT BRAIDING GATES 


As such, the probability of this process not succeeding within 
n attempts is at most 

Prob(xi,...,x n ± 0) < , (37) 


Thus far, we have focused only on measurement and fusion 
operations. We will need to utilize braiding operations as well, 
so we now present the single qubit braiding gates for the 1111 
and 1221 qubit encodings of JK 4 . 
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A. 1111 Single qubit gates 


In the 1111 encoding, we can produce two single qubit 
gates (and any gates generated by these) by braiding adjacent 
quasiparticles. These are given by 

till 



and 

lilt 



till 



E (43) 

c=0,2 


The braiding operation in Eq. (l4ll) gives the computational 
gate 


R-k/3 


1 0 

0 w 


(44) 


where u> = e' ^ . The braiding operation in Eq. < l42b gives the 
gate 




1 

H )\/2 


—LO 


(45) 


The gate set {R^/ 3 , G} is not a computationally universal sin¬ 
gle qubit gate set, since it generates a finite set of 12 gates. 
This set of gates forms a projective representation of the alter¬ 
nating group A4. 


B. 1221 Single qubit gates 


and 

1221 



1221 



E [Fl 22 ] ac Rf [(F? 22 )- 1 ]^ (48) 

c—0,2 


The pure braid operation in Eq. (l46l > gives the computa¬ 
tional gate 


Z = 


1 0 
0 -1 


(49) 


and the braiding operation in Eq. ( l47l > gives the computational 
gate 



1 -iV 3 

-iy/3 1 


(50) 


We note that exchanging the two charge 1 quasiparticles sim¬ 
ilarly generates the gate B. 

The gate set { Z, B} is not a computationally universal sin¬ 
gle qubit gate set. Indeed, it generates a finite set of 6 gates, 
which forms a projective representation of the permutation 
group S 3 (which is equivalent to the dihedral group D 3 ), since 
Z 2 = —B 3 = 1 and BZ = ZB- 1 . 


V. FUSION-BASED NOT GATE 


We now describe a protocol that generates the NOT gate 


X = 


0 1 
1 0 


(51) 


on a 1221 qubit through simple fusion operations that involve 
no measurements. Starting from a qubit in the 1221 encoding, 
we follow the steps: 


In the 1221 encoding, we may also produce unitary single 
qubit gates by braiding. However, in this case, one cannot 
use all braiding exchanges of adjacent quasiparticles, because 
their topological charges are not all identical. We must restrict 
to the braiding transformations which return the quasiparticles 
to the initial configuration of charges. The gates that may be 
obtained from such operations are generated by the following 



1. Introduce an ancillary pair of charge 4 quasiparticles 
that are pair-produced from vacuum. 

2. Fuse one of the charge 4 quasiparticles with one of 
the charge 2 quasiparticles and fuse the other charge 4 
quasiparticle with the other charge 2 quasiparticle. 

This process is described diagrammatically by 
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where -ia = 4xa = 4 — a indicates that the effect of this 
process is the application of a NOT gate. The computational 
gate resulting from this process is independent of the details of 
how the charge 4 quasiparticles are introduced and/or braided 
around the other quasiparticles, because the charge 4 quasi¬ 
particles are Abelian, and so these details can only result in 
unimportant overall phases on the state. Moreover, no mea¬ 
surements are necessary (except, perhaps, to ensure that the 
ancillary quasiparticles indeed carry topological charge 4), be¬ 
cause the fusion outcome involving a quasiparticle of charge 
4 is unique, i.e. 4x2 = 2. We recall that the NOT gate may 
be used to take us between the encoding changing operations 
p(t) to p( 3 ) of Sec. UITBl 

The gate set { X, Z, B} is not a computationally universal 
single qubit gate set, either. Indeed, it generates a finite set 
of 12 gates, which forms a projective representation of the 
dihedral group Dq. 


Steps 1 and 2 constitute a semi-forced (repeat until success) 
measurement, in the following sense. The measurement out¬ 
come z = 2 is undesired, but, as we will explain, it may be 
“undone” in a manner similar to the forced measurement pro¬ 
cedure, since this measurement outcome does not collapse or 
otherwise alter the encoded state. In contrast, the measure¬ 
ment outcomes z = 0 and 4, will collapse the state, as we will 
also see. The measurement outcome z = 4 is not necessarily 
desirable, but it is generally not possible to undo the measure¬ 
ment in this case, since this measurement outcome applies a 
projection to the encoded state. Thus, we must treat measure¬ 
ment outcomes z = 0 or 4 as the “desired” outcomes and end 
the semi-forced measurement process once such an outcome 
is achieved. This process may take n attempts. 

Diagrammatically, the nth iteration of step 1 is given by 


VI. TOPOLOGICAL QUBIT FUSION 


Ref. [Q introduced protocols, referred to as “topological 
qubit fusion” (TQF), which act on multiple topological qubits 
in a manner that reduces the number of topological qubits 
through a series of fusion and measurement operations. The 
reduction of the computational state space using TQF may oc¬ 
cur in a nontrivial way, in the sense that it is not simply a 
projection applied to one of the qubits, but rather may have 
an effect which is equivalent to applying entangling gates to¬ 
gether with measurement projections. TQF was demonstrated 
to be a useful protocol for converting ancillary states into com¬ 
putational gates m, which is what we will use it for in this 
paper. 

We now describe a TQF process for two qubits in the two 
qubit 122221 encoding (see also Ref. |[28i]), reducing a pair of 
qubits to a single 1221 qubit. If we wanted to start from two 
qubits in the 1221 encoding, then the preliminary step would 
be to apply the protocol of Sec. 1111 Cl to switch into the two 
qubit 122221 encoding. 

Starting from two qubits in the two qubit 122221 encod¬ 
ing (enumerating the quasiparticles 1-6 from left to right), we 
follow the steps: 

1. Measure the collective topological charge z € {0, 2,4} 
of quasiparticles 3 and 4. If z — 2, go to step 2. If 
z = 4, go to step 3. If z = 0, go to step 4. 

2. Measure the collective topological charge v € {1,3} of 
quasiparticles 4, 5, and 6. Go to step 1. 

3. Introduce an ancillary pair of charge 4 quasiparticles. 
Fuse one of the ancillary charge 4 quasiparticles with 
quasiparticle 4 and the other one with quasiparticle 5. 

4. Remove the (now ancillary) quasiparticles 3 and 4 in 
a manner that does not create entanglement, e.g. fuse 
them together into vacuum or transport them as a pair 
along the same path away from the other quasiparticles. 



where v\ = 1. The probability of measurement outcome z n = 
2 at the 7/th attempt is always p 2 = independent of the 
encoded state. This follows from the fact that [-F b ° 22 ] 12 = 
and [.F b a22 ] 32 = —for all a, b £ {1, 3}. We emphasize that 
this verifies our claim that this measurement outcome does not 
collapse the state encoded in the two qubit degrees of freedom 
a and b. Thus, the probability of achieving the desired result 
z n = 0 or 4, is po + P 4 = \ at each attempt, though we will 
see that the individual probabilities p 0 and p 4 will depend on 
the encoded state. It is easy to see that the fusion rules require 
that a = b when 2 = 0 and that a = -<b = 4 — b when 
2 = 4 (and that 2 = 2 does not impose any additional relation 
between a and b). 


The nth iteration of step 2 (which is only used when z n = 
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2 ) is diagrammatically given by 


122221 



a ft 


1 2 2 2 2 1 





(54) 


Again, the outcome is independent of the encoded state and 
does not collapse nor alter it. 

Step 3 simply decouples the quasiparticles 3 and 4 from 
the other computational quasiparticles when the measurement 
outcome is z = 4. The process in step 3 (which is only used 
when z n = 4) is represented diagrammatically by 



Step 4 removes the now ancillary quasiparticles 3 and 4, 
leaving a single qubit in the 1221 encoding. 

We denote the operation corresponding to the TQF process 
described in this section, with final measurement outcome z — 
0 or 4, by Q^ z \ with the effect on the two qubit basis states 
given by 

\a,b)^Q^\a,b)= £ Q { :l b \c) (56) 

c—1,3 


where 

«£,* = 22 L [ p r 

Evaluating the F-symbols, we have 

b=c;o,b — -^°a,bxzOac, 

which we can write in matrix notation as 


Q<°> = 


75 0 0 0 
0 0 0 71 


Q (4) = 


0 0 


(57) 

(58) 


0 75 0 0 


V5 


0 


( 59 ) 


Thus, for a two qubit state |T') = b _ 1 3 'T ai b|a, b), the 

effect of the TQF process with measurement outcomes z = 0 
and 4 are, respectively, given by the maps 


l*> ^ 


Q (0) |*> *i,i|1)+®3,3|3) 

vW + fci 2 


Q (4) |ft) = all) + ^3 ,i|3) 

H® (4) l $ >ll ^/l^ 1 , 3 ! 2 + l^aal 2 


(60) 

(61) 


to single qubit states. The possible outcomes z = 0 and 4 of 
the entire TQF process will, respectively, occur with probabil¬ 
ities po = |\ET,i | 2 + I^Q^I 2 andp 4 = |'f'i i3 | 2 + 1 ^ 3 ,i| 2 - We 
emphasize that this TQF process is not deterministic, as the 
probability of these outcomes generally depend on the quan¬ 
tum state encoded in the qubits. As such, they generically 
collapse the initial encoded state (in some nontrivial fashion), 
so one should be careful how and when the TQF process is 
used. 


A. Converting states into gates 


As demonstrated in Ref. IH, one of the useful applica¬ 
tions of TQF is to convert ancillary states into operations on 
a computational state. If we use an ancillary qubit in a state 
with equal magnitude superposition of its basis states, this will 
generate a unitary phase gate on the computational qubit. If 
we use a non-balanced ancillary qubit state, the effect on the 
computational state will generically not be unitary, since it 
will involve some projection. 

We now describe this protocol for JK 4 anyons with qubits 
in the 1221 encoding and the ancillary qubit in the state 

( e | 1 ) + e ^\ 3)) ■ (62) 

Starting from a computational state |'T) with a qubit (at 
least one) in the 1221 encoding and an ancillary qubit state 
|i? 0 / 2 ) in the 1221 encoding, we follow the steps: 

1. Switch the two 1221 qubits (computational and ancil¬ 
lary) into the two qubit 122221 encoding. 


2. Apply TQF to these two qubits. 

This protocol will have TQF outcome z = 0 and 4 with 
equal probability po = P4 = Using the analysis of our TQF 
operation, the resulting operations for these two outcomes are, 
respectively 

|*>| R 4 / 2 ) ^ v / 2Q (0) |'T)|i^/2> = i^/ 2 |*> (63) 

\*)\R H 2 ) ^ v / 2Q (4) |4')|i? 0 / 2 )=i?_ 0/2 |T). (64) 

We emphasize that this protocol essentially consumes the an¬ 
cillary qubit in order to convert it into the phase gate 


R±4>/2 — 


1 0 

0 e ±i(t> 


(65) 


acting on the computational state, where the + and — phases 
are obtained with equal probabilities. 
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VII. IRRATIONAL PHASE GATE 


qubits in the states 


The topologically protected computational gates obtained 
from braiding transformations in the 1221 qubit encoding, 
which are generated by the gate set {Z, B), are not compu¬ 
tationally universal. (Supplementing the braiding gates with 
the gate X does not make them computationally universal, 
either.) We propose to supplement these operations with an 
irrational phase gate 


K — R a /2 — 


1 0 

0 e ia ’ 


( 66 ) 


l$+i,i> 


I < E > —1,3) 



( 68 ) 

(69) 


We generate ancillary qubits in these states using the follow¬ 
ing protocol. 

Starting from an ancillary qubit in the 1111 encoding (enu¬ 
merating the quasiparticle 1-4), we follow the steps: 


where e*“ = . ~ 1 +* 4 V^ - Olmsted’s theorem ll29h assures us 
that a/27r is an irrational number |[451 and, hence, K has infi¬ 
nite order. 

It is straightforward to demonstrate that the gate set 
{ X , Z, 11. K} is computationally universal for single qubits. 
In order to verify that any given set of single qubit gates 
is computationally universal (dense in the set of all single 
qubit gates), we can use the fact that the only infinite proper 
Lie subgroups of PU(2) = SO(3) are isomorphic to either 
U(l) = SO(2) or to 0(2). It follows that, if an infinite sub¬ 
group of PU(2) is not isomorphic to a subgroup of either U(l) 
or 0(2), then its closure must be equal to PU(2). Clearly, the 
gates B and K projectively generate an infinite subgroup of 
PU(2). However, II and K do not (projectively) commute 
with each other, so they cannot generate a subgroup isomor¬ 
phic to U(l). Additionally, we notice that the (nested) group 
commutator [\B, AT -1 ], [B, AT]] ^ 1, where the group com¬ 
mutator [a, b] is defined here to be a -1 fr ] ah. (There is no 
distinction between checking this simply by multiplying ma¬ 
trices and checking it in the projective quotient SO(3), since 
overall phases cancel inside a commutator.) This implies the 
infinite subgroup generated by B and K cannot be a subgroup 
of any 2-stage solvable group, such as 0(2). Thus, the gate 
set {B, K} must generate a subgroup of PU(2) that is dense, 
i.e. has closure equal to PU(2). 

We have shown in Sec. I VI Al how to convert an ancillary 
qubit in the state | R^/?) into a unitary phase gate R± ( j ) / 2 , so 
what remains to be shown is that we can generate the irrational 
state 


I K) 


V2 


1^ ID 


(e |1) + |3» = 



Ra/2/ 


13) 

(67) 


in the 1221 encoding, using measurements, braiding, and fu¬ 
sion operations. The previous sections have developed all the 
operational tools we need to produce this state, so it only re¬ 
mains to assemble them. 


A. Generating the irrational state K) 

As a preliminary step in producing an ancillary 1221 qubit 
in the irrational state | AT), we first prepare two ancillary 1221 


1. Initialize the 1111 qubit in the state |2). 

2. Perform a braiding operation interchanging quasiparti¬ 
cles 2 and 3. 


3. Switch from the 1111 encoding to the 1221 encoding, 
using the protocol of Sec. lIIIBl 

Step 1 may be carried out in various different ways. For 
example, one may initialize the qubit in the state |2) using 
a forced measurement essi. Specifically, this means first 
measure the collective topological charge of quasiparticles 2 
and 3 (the outcome 0 or 2 is unimportant) and then measure 
the collective topological charge of quasiparticles 1 and 2. If 
the outcome of the second measurement is charge 2, then the 
desired state has been prepared. If the outcome of the second 
measurement is charge 0, then simply repeat these two mea¬ 
surements until the outcome of the second measurement is 2. 
Another way to obtain this initialized state is given by replac¬ 
ing the measurements of quasiparticles 2 and 3 in this forced 
measurement procedure with braiding operations of quasipar¬ 
ticles 2 and 3. 

In step 2, the braiding operation may be counterclockwise 
or clockwise, which we will distinguish with a label s = +1 
or —1, respectively. This applies the gate G s to the qubit. 

In step 3, the encoding change involves a measurement out¬ 
come x = 1 or 3, as described in Sec. I III BI If the desired 
measurement outcome x is not obtained, we can always apply 
a NOT gate to switch to the desired outcome, as previously 
explained. There is also a second measurement in this step, 
which is assumed to be y = 2, in order to result in a state in 
the 1221 encoding. When this measurement yields y ^ 2, we 
must discard or recycle the quasiparticles and start over. This 
is not a problem, since we are only generating ancillary states 
at this stage. 

The resulting ancillary qubit produced from this protocol is 
in the 1221 encoding and in the state 


P(*)G S | 2 ) 

s ’*> - ||P(*)GS|2)||' 


(70) 


It is straightforward to check that this gives the states in 
Eqs. (1681) and (l69l ). 

Finally, in order to obtain an ancillary 1221 qubit in the 
irrational state |AT), we start from two ancillary 1221 qubits, 
respectively prepared in the states |$ +1-1 ) and |$_! 3 ), and 
follow the steps: 
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1. Switch the two 1221 qubits into the two qubit 122221 
encoding. 

2. Apply TQF to these two qubits. 

From the analysis of Sec.IVll we find that the result of the TQF 
operation on this pair of ancillary qubits with measurement 
outcome z = 0 is the desired ancillary 1221 qubit state 


|^ > +I,l)|^ > -1,3) H► 


Q (0) |S +1 ,1>I*-1,3> 

||Q(°>|s+i, 1 >I*-i,3>|| 


I K). (71) 


When the TQF operation has measurement outcome z = 4, 
the resulting state is not in the desired form, nor can it easily 
be salvaged for our purposes. Since these are ancillary qubits, 
we may simply discard or recycle the final qubit when the 
TQF operation has z = 4. 


B. Converting | K) into the phase gate K 

Once we have an ancillary qubit in the 1221 encoding pre¬ 
pared in the irrational state | AT), we can use TQF to convert 
the ancillary qubit into a phase gate acting on a computational 
qubit in the 1221 encoding, following the protocol described 
in Sec. IVI Al However, this is not a deterministic process, as 
the measurement outcomes z = 0 and 4 of the TQF procedure 
will occur with equal probabilities Po = P 4 = The out¬ 
come z = 0 results in an application of the unitary phase gate 
AT to the computational state, while the outcome z = 4 results 
in an application of AT -1 . If we intended to apply a K gate 
in our quantum computation, but instead generated the gate 
A' -1 from this protocol, we need a strategy for correcting this 
undesired outcome. We obviously cannot correct this simply 
by applying gates from our deterministic gate set {X, Z, 11). 

A simple strategy for dealing with this issue is to repeat¬ 
edly apply the protocol for converting states | AT) into phase 
gates until the product of gates applied is equal to the desired 
phase gate K. More specifically, after one application of the 
protocol, we will generate the AT gate with probability 1/2 
and A' -1 with probability 1/2. If we generated the undesired 
gate AT -1 , then we apply the protocol two more times. These 
two applications have probability 1/4 of generating K 2 , prob¬ 
ability 1/2 of generating 1, and probability 1/4 of generating 
AT -2 . Thus, following the initial gate AT -1 , the two additional 
applications of the protocol have probability 1/4 of making 
the total product of gates generated equal to K. If the desired 
outcome is not achieved after these two additional applica¬ 
tions of the protocol, then we continue to repeatedly apply 
the protocol in this way, until the desired outcome is finally 
achieved. 

We can think of this as carrying out a random walk on the 
integers, where the position £ Z after n steps is the expo¬ 
nent of the product K Xn of gates after n applications of the 
protocol. As such, the random walk starts from xo = 0 and 
each step has an equal probability of moving the position to 
x n +i = x n + 1 or x n - 1. 

This random walk strategy for generating the K gate by 
repeatedly converting | AT) states into phase gates using TQF 


until successful is only viable if the probability of achieving 
the desired result goes to 1 as the number n of repeated ap¬ 
plications of the protocol gets large (n —> oo). We analyze 
the random walk in Appendix [B] and find that this is indeed 
the case, but that the probability of not achieving the desired 
result (a positive valued position) in n steps, where n is odd, 
is 

fl \j 

Prob (xi, ■ ■ ■, x n < 0) = - —(72) 

(n + 1)!! 


which goes to zero as y ^ for n large. This is certainly 
less ideal than an exponentially fast convergence (which is 
what we found for all our other probabilistic protocols), but 
it nonetheless permits a higher level strategy that allows the 
use of such random walk generated gates to be used in BQP 
(bounded error quantum polynomial time) quantum computa¬ 
tions. 

A strategy for using such random walk generated AT gates 
while satisfying the conditions for BQP is the following. If 
the quantum computation we wish to carry out involves k ap¬ 
plications of the gate AT, then we allot each K gate at most 
k 2 steps in the random walk used to generate it. For each K 
gate, if the random walk reaches k 2 steps without achieving 
the desired result, we terminate the process and consider the 
gate and, hence, the entire computation to have failed. Each AT 
gate that we attempt to implement in this way will thus have a 
probability 


PK -fail 


(fc 2 )!! 

{k 2 + 1)!! 


(73) 


of failing and 1 — PK-m of being successfully generated. 
Consequently, the entire computation will have a probability 
(1 — PK-M\) k of being successfully implemented (at least with 
respect to the AT gates’ random walk issue). In the limit as k 
gets large, we see that the probability for the entire computa¬ 
tion being successfully implemented is 


lim (l-p K -fmi) k = lim 

k—t oo k—>oo 




Thus, the probability of failure of the computation due to the 
random walk generated K gates using this implementation 
strategy is bounded, and hence this strategy satisfies the con¬ 
ditions for BQP. (We note that it can be inferred from the mea¬ 
surement outcomes whether or not the AT gate has been suc¬ 
cessfully applied, so there is no requirement for the success 
probability to be greater than 1/2.) However, we notice that 
the number of operations that may be required to implement 
all k applications of K gates using this strategy scales as k 3 , 
so the polynomial exponent of the computation length scaling 
increases by a factor of 3. 

Clearly, the crude strategy described here is far from ideal, 
as our goal was merely to demonstrate the existence of a strat¬ 
egy that works, in principle, within BQP. There are a num¬ 
ber of strategies for algorithm optimization. One simple way 
to reduce the negative impact (i.e. the costliness and conver¬ 
gence issues) of the irrational AT gate is to compile the quan¬ 
tum algorithm with an aversion to the K gate. In other words, 
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when synthesizing gates or algorithms from the generating 
gate set {X, Z, B 1 K}, the optimization should minimize a 
weighted combination of both the word length and the total 
number of K gates used. 

Another strategy for improving the situation is to employ a 
sliding cutoff for the number of attempts allowed before termi¬ 
nating the random walk attempting to implement a K gate. In 
particular, one may wish to terminate random walks sooner at 
earlier stages of the computation, while allowing longer walks 
as one nears the end of the computation. A more drastic modi¬ 
fication of the strategy is to pick a cutoff point in each random 
walk (which could even be after one step), after which one re¬ 
compiles the algorithm to see if a different path forward would 
be more economical 0. 


VIII. CONTROLLED-Z ENTANGLING GATE 


A. Generating entanglement resources 

We wish to generate ancillary pairs of qubits in the 1221 
encoding which are in the entangled two qubit state 

\*h) = H\$ + ) = \ (111) + 113) + |31) - 133)). (77) 

[We use the conventional notation | < E >± ) = ^(|11) ± |33)) 
and |4' ± ) = ^=(|13) ± |31>) for the Bell states.] The 
Hadamard operator here can be equivalently applied to either 
the first or the second qubit in this expression. We now pro¬ 
vide a protocol for generating such states. 

Starting from two ancillary qubits in the 1221 encoding 
(enumerating the quasiparticles 1-8 from left to right), we fol¬ 
low the steps: 


The set of single qubit gates { X, Z, B, K } is computation¬ 
ally universal for single qubits. If we can supplement this gate 
set with any entangling gate, we would obtain a full computa¬ 
tionally universal gate set. In this section, we provide a proto¬ 
col for generating the two qubit entangling gate, controlled-Z 


C(Z) 


■ 1 0 0 0 - 

0 10 0 

0 0 10 

.000 - 1 . 


(75) 


In addition to the inexpensive gates and operations (such 
as topological charge measurements, braiding, and TQF), our 
protocol will utilize a single application of the Hadamard gate 


V2 



(76) 


1. Initialize each 1221 qubit in the state |+) = H\l) = 

*(|1> +13))- 

2. Measure the collective topological charge r £ {0, 2,4} 
of quasiparticles 3, 4, 5, and 6. If r ^ 0, go to step 1. If 
r = 0, go to step 3. 

3. Perform a braiding operation that interchanges quasi¬ 
particles 3 and 4 and one that interchanges quasiparti¬ 
cles 5 and 6, with the opposite chirality. 

4. Move the quasiparticles 3, 4, 5, and 6 to the right of 
quasiparticles 7 and 8 in a manner that does not create 
entanglement (i.e. have them all follow the same path). 

5. Apply the Hadamard gate H to one of the qubits. 


The Hadamard gate is costly, as it requires a long string of 
the generating gates to accurately produce, including a large 
number of applications of the irrational phase gate K, each 
of which is relatively expensive to implement, as seen in 
the previous section. More specifically. Ref. [3l|] developed 
algorithms for synthesizing single qubit gates from the set 
{Z, B, K} that produces approximate implementations of H 
to precision e with /\-count around 41og 7 (l/e). For exam¬ 
ple, they produce an approximate H gate with trace preci¬ 
sion e < 10~ 8 using a string of 67 gates, 40 of which are K 
gates 0. 

However, the single application of H is utilized during the 
ancillary state preparation stage of the protocol, so its gener¬ 
ation can be carried out in a parallelized fashion that avoids 
the previously discussed issues with applying K to the com¬ 
putational state within a quantum computation. The use of the 
Hadamard gate also means that the precision of our C(Z) gate 
will depend upon the precision to which we approximate the 
Hadamard gate H that is applied in this protocol. 

We split the protocol for generating C(Z) into two parts: 
(A) generating the ancillary two qubit states that will serve as 
entanglement resources, and (B) converting the entanglement 
resource ancillary state into the C(Z) gate using TQF. 


Step 1 may be carried out by applying a Hadamard gate H 
to the qubit initialized in the basis state |1). However, this 
is not the best way to perform this initialization, since the 
Hadamard gate is costly and not exact. A better method is 
suggested by Eq. (1231) and the corresponding F-symbol 




ab 


l 

71 


i i 
i -i 


= H ab 

ab 


(78) 


where a = 1, 3 and b = 0, 2. From this, we see that, for a 
qubit in the 1221 encoding, we can simply measure the topo¬ 
logical charge b £ {0, 2} of the pair of charge 2 quasiparti¬ 
cles. (We could, equivalently, measure the topological charge 
b £ {0, 2} of the pair of charge 1 quasiparticles, if desired.) If 
the measurement outcome is b = 0, then the 1221 qubit is in 
the desired state |+). If the measurement outcome is b = 2, 
then the 1221 qubit is in the state |—) = Z |+), so we merely 
apply the gate Z to obtain the desired initial state. Similarly, 
we could initialize the 1221 qubit in the state |+) by pair pro¬ 
ducing a pair of charge 1 quasiparticles from vacuum and pair 
producing a pair of charge 2 quasiparticles from vacuum and 
then aligning them into the 1221 qubit configuration. 









14 


The measurement in Step 2 is diagrammatically given by 


12211221 



a b 


12211221 



a b 


When r = 0, this is equal to 


12211221 



Steps 3 and 4 are simply rearranging the configuration of 
quasiparticles to that of two qubits in the 1221 encoding, in a 
manner that does not affect the state. Now, the original quasi¬ 
particles 1, 2, 7, and 8 comprise the first 1221 qubit and quasi¬ 
particles 4, 3, 6 , and 5 comprise the second 1221 qubit. 

Combining steps 2, 3, and 4 (with measurement outcome 
r = 0 ), the effect on basis states can be written diagrammati¬ 
cally as 


12211221 



a b 

12211221 



The effect on the initialized states (i.e. combining steps 1-4) 
is thus 

|+}|+) ^ |$+). (82) 

Step 5 takes us from the Bell state |<1> + ), obtained by apply¬ 
ing steps 1-4, to the desired state 

|$ ff ) = 1 (g> iT|$+) = H® 1|$+). (83) 

In contrast with the initialization step 1, we do not know a way 
to circumvent the use of the costly Hadamard gate H in this 
last step. 

B. Converting \& h ) into a controlled-Z gate 

Assuming that we now have two ancillary qubits in the 
1221 encoding paired up in the entangled two qubit state 


<!>//), we can convert this ancillary state into the C(Z) gate 
acting on two computational qubits in the 1221 encoding us¬ 
ing the following protocol. 

Starting from two computational qubits (A and B ) in the 
1221 encodings and an ancillary pair of qubits (1 and 2 ) in 
the 1221 encoding in the entangled state < 1 >//), we follow the 
steps: 

1. Perform TQF on qubits A and 1, with measurement out¬ 
come zai £ {0,4}. If zai = 4, apply the gate Z to 
qubit A. 

2. Perform TQF on qubits B and 2, with measurement out¬ 
come Zb2 £ {0,4}. If zb2 = 4, apply the gate Z to 
qubit B. 

It is straightforward to check that the effect of this protocol, 
which can be written for the computational qubit basis states 
as 

\a,b)\* H ) ^ 2Q i 2lQ ( §l\a,b)\^ H ) (84) 

\a,b)\<f> H ) IZbQ^bIM^h) (85) 

\a,b)\*H) ^ 2Z A Q < £ 1 Q ( gl\a,b)\$ H ) ( 86 ) 

M)|$h> 2Z A Z B Q%lQ%l\a,b)\9> H ) (87) 

for the four possible combinations of TQF measurement out¬ 
comes zai and zb 2 , all give the same result, which is an ap¬ 
plication of the controlled-Z gate to the two computational 
qubits 

\a,b)\$ H )^C(Z)\a,b)- (88) 

We emphasize that this protocol consumes the ancillary entan¬ 
glement resource pair of qubits in order to convert them into 
the C(Z) gate acting on qubits A and B of the computational 
state. 


IX. CONCLUSION 

In the example of SU( 2)4 and JK 4 anyons, we have seen 
that it is possible for an anyon model to not be computa¬ 
tionally universal with braiding alone, but to become compu¬ 
tationally universal when braiding is supplemented with fu¬ 
sion and measurement operations. This demonstration of the 
utility of fusion and measurement operations encourages fur¬ 
ther analysis of such operations and exploration for additional 
theories whose computational power can be supplemented or 
even made computationally universal in this way, particularly 
for theories that would be easier to physically realize. 

The same strategy of supplementing braiding operations by 
fusion and measurement operations can also be considered 
for theories with symmetry defects, which are described by 
G-crossed braided categories, in which braiding is general¬ 
ized to incorporate symmetry action [ij]. The defect theories 
that are relatively easy to physically realize appear to have 
(G-crossed) braiding which is not computationally universal, 
so it would be interesting to determine whether any of them 
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could be made computationally universal via fusion and mea¬ 
surement operations. One known example of a defect theory 
which may benefit from fusion and measurement operations is 
a bilayer Ising TQFT system with layer interchange symmetry 
(■G = Z 2 ); fusion and measurement may be used to switch be¬ 
tween a qubit state being encoded in a quasiparticles, where 
they may be relatively easier to physically manipulate, and 
being encoded in defects (also known as “genons”), which 
provide computational universality through braiding |33| . 

Appendix A: Basic Data of the Anyon Models 

In this section, we provide general expressions for the basic 
data (topological charges, fusion rules, /-’-symbols, and R- 
symbols) of the SU(2)fc and .1K/, anyon models, for general 
k. For JK 4 , we also evaluate and tabulate the /-’-symbols and 


//.-symbols that we use in this paper, for convenience. 


1. SU (2) fc Anyons 

The SU(2) fc anyon models (for k an integer) are “q- 
deformed” versions of the usual SU(2) representation the¬ 
ory 0 . Roughly speaking, this means integers n are re- 

n/2 _ —n/2 

placed by “g-integers” [n} q = q 1 / 2 _ q _ 1 / 2 » where the defor- 

• 2 7T 

mation parameter q = e l k + 2 is taken to be a simple root of 
unity. These anyon models describe SU(2) fc Chern-Simons 
theories ll35ll and WZW CFTs ll36[ l37ll . and give rise to the 
Jones polynomials of knot theory 113811 . Their braiding statis¬ 
tics are known @1 to be computationally universal for all k, 
except k = 1, 2, and 4. The anyon models are described by: 



where { } ? is a “g-deformed” version of the usual SU(2) 6j-symbols. The sum over z in this expression is over all integers for 

which each term is well-defined, i.e. < z < 2 max , where 2 min = max{ji +j 2 +j 3 +j, ji +j r2 +j 3 + j 23 , ji +J12 + j_ + J23} 
and 2 max = min{ji + j 2 + jv2-,JV2 + j 3 + j,h + jz + J23, ji + J23 +/}■ The Frobenius-Schur indicator = dj 
gauge invariant quantity which plays a role in bending and straightening lines. 


pjjj 

3 


J 00 


2. JKfc Anyons 


There are anyon models based on the Jones-Kauffmanbracket [39,4dl- which may be derived from Temperley-Lieb recoupling 
theory. These anyon models, which we denote as JK/ r , are closely related to the SU(2) fc anyon models. They have the same 
number of topological charges and the same fusion rules as SU(2) fc at the corresponding level, and the basic data shares a 
very similar structure. It is conventional to label the topological charges of JK/, by integers, rather than both integers and half¬ 
integers, so there will be a translation between these by multiplying/dividing the charge labels by 2. We define A = ie ~ l ' 2( - k + 2 '>, 
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and [n\ A = 


A 2n_ A -2n 

A*-A~* ' 


The any on model data is given by: 



The sum over z in this expression is over all integers for which 
each term is well-defined, i.e. z mm < z < z m ax , where 
Zmin = max{a + b + c + d, a + e + c + f, b + e + d + f} 
and z max = min{a + b + e, e + c + d,b + c + f, a + f + d}. 
It is straightforward to check that the conditions for complete 
isotopy invariance are satisfied for JK/,, i.e. all the bending 
transformations in this theory (for the choice of gauge used 

above) are trivial, since [T^* ab ] 0c = [F“' bb ] cQ = for 

all a, 6, c with N£ b ^ 0. We note that k = 2 is the Ising anyon 
model. 


Comparing the basic data of JKfc to that of SU(2)fc, we no¬ 
tice that the structure of these expressions and their SU(2)fc 
counterparts are nearly identical. The main difference is that 
the (/-deformation factor is changed by a minus sign (note 
that A 2 plays the role of g), which requires minus sign fac¬ 
tors that necessarily arise in the associativity expressions, as 
well as factors of i that arise in the braiding. Notice that the 
Frobenius-Schur indicator >c a = d a [F “ oo ] 00 h ere I s trivial; as 
required for the theory to be completely isotopy invariant. 

In fact, the two theories can be related by gluing a semion 
onto the odd/half-integer charges. More specifically, this 
means we take the direct product of one of these theories with 
a semion and then restrict the topological charge set so that 
odd charges in the JK/, sector (or half-integer charges in the 
SU( 2 )/, : sector) and the nontrivial charge in the semion sec¬ 
tor are always paired up. We can write this as SU (2) fc = 

7 ( 3 / 2 ) 


JK fc x Z^ 3/2) 


, where 7 ; = SU (2 ) 1 is the (Z 2 ) semion 


with Fl 1 11 = —1 and 


1 = —i, and C = {(a, b)\a € 
{0,1, 6 £ {0,1}, a = b mod 2}. This relation 

is confirmed by the corresponding values of the topological 
twists and S'-matrices of the respective theories. 


3. Some Useful / -symbols and /7-symbols of JK4 


In this section, we tabulate the ^-symbols and /7-symbols 
of JK 4 that are used for calculations in this paper, for conve¬ 
nience. 



Appendix B: Random Walk 

In this section, we consider a random walk over the inte¬ 
gers Z, starting at position x = 0, with each step having equal 
probability i of taking a step from x to x + 1 or x — I. We 
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compute the probability of never taking a positive valued po¬ 
sition during a walk of n steps or, equivalently, of never taking 
a negative valued position during a walk of n steps. Since a 
walk starting at 0 can only become positive on an odd step, it 
is clear that this probability is the same for n and n + 1 steps, 
when n is odd. As the problem is symmetric between avoiding 
positive and negative positions, we will simplify notation by 
computing the probability of the position never going negative 
within n steps. We compute this iteratively by counting the to¬ 
tal number of paths one can take without ever going negative 
within n steps for n = 2m — 1 odd. 

First we note that the total number of unconstrained paths 
possible after n = 2m — 1 steps is = 2” = 2 2m_1 . 

After the first step (m = 1), there is equal probability of 
going to x = +1 and —1, so the number of paths that end up 
at positive position x without ever going negative is simply 

NW=5x, i. (Bl) 

Since each step has equal probability of moving the position 
by +1 or — 1, a consecutive pair of steps can move the position 
by +2, 0, or —2, with there being one way to move ±2 and 
two ways to remain in the same position (i.e. first step forward 
and then step backward, or first step backward and then step 
forward). Thus, we can write the iterative expression for the 
number Nx™' 1 of possible n = 2m — 1 step paths that end up 
at positive position x without ever going negative as 


of never going negative within n = 2m— 1 steps can similarly 
be expressed iteratively as 


(m+l) 

p+ 


N\ 


1+1) 


4ivj m) - N[ m) 


n: 


(m+l) 


total 


= i- 


N, 


(m) 


AN 


V 

where we used the property 


AN 


(m) 


total 


(m) 


(m) 

P+ 


2 (m + 1) 


(m) 

P+ , 


(B6) 


N[ m) = 2 


-w| m) . 

m+l + 


(B7) 


Solving the resulting iterative expression of Eq. (1B6I) . we ob¬ 
tain the probability of avoiding going negative in n = 2m — 1 
steps to be 


(m) 

P+ 


n 

p=i 



(2m — 1)!! 
(2m)!! 


(B8) 


( jy(m+l) = 2 N (jn.) + N (m) for x = l 

[ Ni m+1) = Njrt + 2 ivi m) + N™ for x > 1 ' 

(B2) 

(We notice that the contribution from position x — 2 is missing 
from x = 1, because that would have included a path that went 
negative.) 

The total number of possible n = 2m — 1 step paths that 
end up at positive position without ever going negative is thus 

=Y J N^ m) , (B3) 

£>0 

which satisfies the iterative expression 

N^ +l) =4N^ n) -N[ m) . (B4) 


In order to obtain the scaling of //"' 1 for m large, we use the 
relations of double factorials to the Gamma function together 
with Stirling’s formula 


(m) 

P+ 


T(m + |) 
v^T(m + 1) 
1 

7^ 


exp mln(m-|—) — (m-|—) + . 


— (m + i) ln(m + 1) + (m + 1) 

1 -I _2s 

m 2 + 0{m 2 3 4 ). 




(B9) 


Thus, the probability 


P+ 


]\r( rr ' 

(m) _ ^+ 


N t 


(m) 


total 


Thus, we find that the probability of taking n steps without 
ever taking a negative valued position goes to zero as n~ 2 as 
(B5) n —> 00 . 
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